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] We analyze the opportunity to observe thermal dileptons emitted off 

deconfined matter resulting in heavy-ion collisions at RHIC and LHC en- 
ergies. Special kinematical conditions provided by the detector systems 
PHENIX and ALICE, and the so-called M± scaling behavior of thermal 
Q\ ' dilepton spectra are taken into account. Our considerations include en- 

<n : 

ergy loss effects of the fast heavy quarks in deconfined matter, which for 
J> ! themselves can help to identify the creation of a hot and dense parton 

in . 

medium. Due to a threshold like effect for decay leptons we find a window 
at large transverse pair momentum and fixed transverse mass where the 

o . 

thermal signal can exceed the background of dileptons from correlated 

Q\ ; 

"*^h , semileptonic decays of charm and bottom mesons. 
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I. INTRODUCTION 



One of the ultimate goals in investigating central heavy-ion collisions at very high 
energies is to analyze the properties of highly excited and deconfmed matter. Direct 
probes, like electromagnetic signals Q , are considered as useful carriers of nearly undis- 
turbed information from the transient hot reaction stages. Pairs of electrons and muons 
are experimentally accessible, and will be measured also in forthcoming heavy-ion ex- 
periments pm . These types of probes have proven useful at Super Proton Synchrotron 
(SPS) energies at CERN in the resonance region below the u, p (§], and in the J/ip 
region ||; also in the so-called continuum region in between a yet unexplained ex- 
cess seems to be found |J. However, recent estimates J7| of dileptons stemming from 
semileptonic decays of open charm and bottom mesons produced in the same central 
collisions of heavy ions show that these represent the dominating dilepton source at 
energies envisaged at the Relativistic Heavy Ion Collider (RHIC) in Brookhaven and at 
the Large Hadron Collider (LHC) in CERN. In particular, the dilepton signal from ther- 
malized, strongly interacting matter in the intermediate reaction stages is estimated to 
be up to two orders of magnitude below the background of correlated decays of open 
charm or bottom meson pairs in a wide invariant mass region. Even strong energy 
loss effects of the heavy quarks in the hot and dense medium [jsHlOU do not basically 
change this situation: The decay dileptons might be suppressed down to the Drell-Yan 
background, but it is questionable whether such strong energy loss effects can happen 



in reality ||11|| . In this situation at least two issues can be seen. First, it is necessary to 
extend the traditional strategy for searching a thermal dilepton signal from deconfmed 
matter. The usual approach relies on the analysis of the invariant mass (M) spectra. 
One can expect, however, that the double differential dilepton spectrum as a function 



of the transverse pair momentum Q± and transverse pair mass M± = J M 2 + Q\ con- 
tains much more information which can be used to impose special kinematical gates for 
finding a window to observe the thermal dileptons in the continuum region in spite of 
the large decay background. Second, one can try to employ the initial hard production 
of open charm and bottom followed by semileptonic decays to probe deconfmed matter 
by energy loss effects. A similar idea for hard jet production is widely under discussion 
now 



12 



In the present paper we focus on both of the above issues. We compare the so-called 
Mi scaling property of the dielectrons stemming from deconfmed matter with pairs 
resulting from correlated open charm and bottom decays. As well-known, under certain 



conditions the dilepton yield from a thermalized quark-gluon plasma depends only on 
the transverse mass, i.e., there is no dependence on the variable Q± at fixed M±_ in 
a given rapidity interval Y. At the same time the corresponding Q± dependence of 
dileptons from heavy quark decays has qualitatively a different behavior: It is peaked at 
small values of Q± and drops rapidly at increasing values of the transverse momentum. 
To get realistic estimates of the M±_ scaling we perform Monte Carlo simulations for 
both the correlated heavy quark decays and the thermal production including the 
acceptance of the detectors PHENIX at RHIC and ALICE at LHC. To strengthen 
the Mi scaling violation for dileptons from decays we select only such leptons with 
individual transverse momentum p± > p™ m with appropriate value of p™ m . We find 
that, while the thermal dilepton yield from the plasma is less affected by the kinematical 
restriction, the contributions from decays at large transverse pair momentum is reduced 
drastically. Actually this leads to a threshold like behavior of the transverse momentum 
spectrum at fixed value of M± for dileptons from heavy quark decays in contrast to the 
approximate plateau of the thermal yield which extends up to the kinematical limit at 
Q± = M±. As a result, an opportunity is found which allows to separate the thermal 
dilepton signal from the background of heavy quark decays on the basis of the M± 
scaling property. We also study to what extent the single electron and dilepton spectra 
are modified by the energy loss effects of heavy quarks when traversing deconfined 
matter. In doing so we concentrate on calculable processes in the first few fm/c of 
the collision course. There are, of course, other processes in the later reaction stages 



that might considerably affect the momenta of the charm and bottom mesons |13| , and 
also thermal radiation of semihard probes can stem from later stages. For instance, in 
ref. [|14| the hypothesis has been formulated, resting on an analysis of present Pb + Pb 
data at SPS energy, that all hadron spectra can be parametrized by a common flow 
velocity and freeze-out temperature. This, if correct for higher energies too, would 
imply that charm and bottom spectra are much stronger modified than assumed in 
the following, and as a consequence the chances to identify a thermal dilepton source 
would be better. 

Our paper is organized as follows. In section 2 we present our calculation pro- 
cedure of charm and bottom spectra and the resulting single electron and dielectron 
spectra. Section 3 is devoted to the thermal dileptons at midrapidity and the M± 
scaling properties. Our conclusions can be found in section 4. 



II. DILEPTONS FROM OPEN CHARM AND BOTTOM DECAYS 



A. Initial charm and bottom production 

In our approach we utilize the leading order QCD processes gg — > QQ and qq — > QQ 
for heavy quark - antiquark (QQ) production and simulate higher order corrections by 
an appropriate Kq factor. We have checked that such a procedure reproduces within 
the needed accuracy the more involved next-to-leading order calculations [|7|. The 
number of QQ pairs, produced initially with transverse momenta p_n = —p± 2 — at 
rapidities yi )2 in central AA collisions can be calculated by 

dN QQ = T A a(0) Kq H(yi, y 2 ,p±) dp\ dy x dy 2 , (1) 
H(yi,y 2 ,p±) = x 1 x 2 <f g (xi,Q )f g (x 2 ,Q 



dt 

+ E Q 2 ) ffa> Q 2 ) + Q 2 ) Q 2 )] d ^t \ , (2) 

qq 



where <j® g /di are elementary cross sections (see for details f7[,|l5|1), fi(x,Q 2 ) with i = 
g,q,q denote the parton structure functions, x\ 2 = m±_ (exp{±yi} + exp{±y 2 }) /y/s 



and mj_ = yp± + ttvq. As heavy quark masses we take m c = 1.5 GeV and m& = 
4.5 GeV. We employ the HERA supported structure function set MRS D'- from 
the PDFLIB at CERN. Nuclear shadowing effects are not included since they will be 
considered separately. In addition, for our present goals we expect no significant mod- 
ification by shadowing in the large transverse momentum region according to results 
of ref. ||17|| . The overlap function for central collisions is Taa{0) — A 2 /(7[R A ) with 
R A = 1.1 A 1 / 3 fm and A = 200. From a comparison with results of ref. we find the 
scale Q 2 = ^Mq and Kq = 2 as most appropriate. 

B. Energy loss effects 

Energy loss effects (cf. |TT]JT^-^TJ and further references therein) have been found 



important p|-|T0| since the momenta of back-to-back moving heavy quarks are degraded 
and the resulting dileptons in correlated semileptonic decays get less invariant mass. 
This process causes a considerable reduction of the number of high invariant mass 
dileptons from correlated DD and BB decays. 

To model the energy loss effects of heavy quarks in expanding matter we assume that 
a heavy quark has the same rapidity as the longitudinal flow at the space-time point of 
the initial quark production. Therefore, with respect to the fluid's local rest frame the 



heavy quark has essentially only a transverse momentum p± which may depend on the 
proper local time r in accordance with the energy loss dE/dx in transverse direction. 
This picture implies a very fast thermalization of gluons and light quarks. 

Unfortunately a sufficiently general treatment of the energy loss rate dE/dx for 
heavy quarks passing through expanding QCD matter has not been developed so far. 
Recently the total energy loss of a high-energy parton propagating transversely through 
expanding deconfmed matter has been derived as dE / dx\ expiin ding = £dE/dx\x f [ PI , 
supposed the medium cools according to a power law and the initial time is short. 
The numerical factor £ is 0(2) for a parton created inside the medium. Within such 
a scenario the stopping power depends only on the final temperature at which the 
deconfmed medium is left, as indicated by the subscript Tf. We denote this stopping 
scenario as model I and take dE/dx\T f from |TT|JT9| --pT| . This model is appropriate for 
LHC conditions, where the initial temperature is very high and where the life time of 
deconfmed matter is large enough that most energetic heavy quarks can leave this stage 



before T c is reached. For RHIC conditions we adopt the results of ref. |19j for dE/dx in 
a static medium but apply them in the simplest generalization to a medium with time 
dependent density. This gives evolution equations for the transverse momentum p±(r) 
of quarks propagating the distance r±(r) in the transverse direction. The scenario is 
denoted as model II, which is described in more detail in refs. 0,0- Strictly speaking, 
due to the finite size of the medium, the large mass of the heavy quarks and other 
approximations inherent in such approaches, the models I and II provide at best semi- 
quantitative estimates. Nevertheless one can hope to cover various limiting cases to 
get some insight in the importance of the energy loss effects. 

To get the spectra of charm and bottom quarks after energy loss we use a Monte 
Carlo simulation with a uniform distribution of the random initial position and ran- 
dom orientation of QQ pairs in the transverse plane. In model I the total energy loss 
is determined by the transverse distance d of a created heavy quark to the boundary 
of the system, the quark initial energy and the temperature when leaving the system. 
Most quarks experience an energy loss oc d and only a few ones ocd 2 . (The transverse 
expansion at early times can be neglected, and we take as averaged transverse radius 
of parton matter R = 7 fm.) In model II we integrate the evolution equations for the 
heavy quark momentum p±(r) and position r±(r) together with the time evolution of 
the temperature T(t) and fugacities A g 9 (r) (see below) as described in refs. @,fnj. The 
quantity rj_(r) is used to check whether the considered heavy quark propagates still 
within deconfmed matter; if it leaves the deconfmed medium it does not longer experi- 



ence an energy loss. (The energy loss in a mixed and hadron phase might be incorpo- 



rated in line with the approach of We consider a heavy quark to be thermalized 
if its transverse mass m± during energy loss becomes less than the averaged transverse 
mass of thermalized light partons at given temperature. These heavy quarks are redis- 
tributed according to the transverse mass distribution dNq/dm^ oc exp{— mj_/T e s} 
with an effective temperature T eS = 150 MeV [ITU . 



C. Single electron spectra 

We employ a delta function like fragmentation scheme for the heavy quark con- 
version into D and B mesons. It results in the same transverse momentum of the 
heavy meson as the heavy parent quark has had previously. Inclusive single transverse 
momentum spectra of electrons/positrons from semileptonic open charm and bottom 
decays, i.e., D(B) — > e + X, and dilepton spectra from correlated decays as well, i.e., 
D{B)D{B) — > e + Xe~X, are obtained from a Monte Carlo code which utilizes the in- 
clusive primary electron energy distribution as delivered by JETSET 7.4. The heavy 
mesons are randomly decayed in their rest system and the resulting electrons then 
boosted appropriately. The average branching ratio of D(B) — > e + X is taken as 12 
(10)% p2[ . The neglect of secondary electrons is justified since we consider here single 
electron spectra with p± > 1 GeV and dilepton spectra with M > 2 GeV. The average 
energies of the secondary electrons, in particular from B decays, are too small to affect 
the spectra in these kinematical regions noticeable. 

The following acceptance cuts are utilized in our Monte Carlo calculations. 
PHENIX [[| can register electrons in the pseudorapidity window \r]\ < 0.35 with az- 
imuthal coverage 0-90 and 135 - 225 degrees, while ALICE || is azimuthally symmetric 
and covers \r]\ < 0.9. As minimum transverse electron momentum we use p± = 1 GeV. 
A maximum transverse momentum cut is not imposed, even at ALICE there is prob- 
ably one. We consider here only the electron and e + e~ channels and the midrapidity 
region since the energy loss effects are strongest there [TTJ . 



The results of our calculations of the transverse momentum spectra of inclusive 
single electrons at RHIC conditions are depicted in figs, la and b. As already pointed 



out in ref. |]2"3" |, a measurement of single electrons with p± > 1.5 GeV at RHIC can give 
a clean and background-free charm signal in heavy-ion collisions. Energy loss effects 
of the heavy quarks modify this signal. As seen in figs, la and b such modifications, 
according to the energy loss described above, are measurable. For both the charm 



and the bottom decay the transverse single electron spectra become softer. This leads 
to a noticeable change of the slope of the spectra. One can fit the distributions by 
dN e -/dp± oc exp{— p±/T e } in the interval 2.5 GeV < p± < 4.5 GeV and finds T e = 
646 (512) MeV in case of no stopping (stopping according to model II) for charm, 
and T e = 851 (640) MeV for bottom. The superposition of charm and bottom decay 
spectra can be parametrized by T e = 711 (574) MeV. Therefore, in comparing the 
p± spectra of single electrons in pp and AA collisions at the same beam energies one 
can reveal the stopping effect. (A direct measurement of the D(B) p± distributions 
would be more useful, but it is unclear whether the D(B) mesons can be identified via 
their hadronic decay modes in such high-multiplicity environments.) In extreme case, 
when the stopping is so large that thermalization with the surrounding matter happens, 
then a pronounced "thermal bump" at p± < 2.5 GeV appears. All of these considerable 
modifications of the transverse electron spectra in AA collisions in comparison with pp 
collisions can be helpful for identifying the creation of a hot and dense parton system. 
It should be noted that the energy loss effects affect the single-electron pj_ spectra 
much stronger than another fragmentation scheme (e.g. the below described Peterson 
scheme) would do. 



D. Decay dielectrons 

The invariant mass dilepton spectra resulting from both correlated charm and bot- 
tom decays are displayed in figs. 2a and b for RHIC and LHC energies. Also shown 
is the distribution of thermal dileptons according to eq. (9) below. As pointed out 
earlier the M spectra are by far dominated by the dileptons from correlated open 
charm and bottom decays. A strong energy loss according to the above models reduces 
the decay contributions, so that (i) at RHIC charm and bottom decay dileptons are in 
the same order of magnitude as the Drell-Yan background |p|,[10[] , and (ii) at LHC still 
bottom will dominate the Drell-Yan background, and charm is slightly below bottom 
0. While without any detector acceptances usually the predictions of the thermal 
yield (dealt with in the next section) is in the same order of magnitude as the Drell- 
Yan background [|7^ 15], 21 1 , it turns out that both the acceptances of PHENIX and 



ALICE suppress the thermal signal stronger than the Drell-Yan background. At the 
first sight such a situation looks quite unfavorable for a measurement of the thermal 
signal even at such high energies as achieved at RHIC and LHC. Nevertheless we show 
below that one can select special kinematical conditions superimposed to the detector 



cuts to analyze the more informative double differential dilepton spectrum with the 
aim to find a window for the thermal signal. 



III. THERMAL DILEPTON PRODUCTION AND M ± SCALING 

A. Spectra of thermal dileptons 

In order to estimate the thermal dilepton yield from deconfined matter we restrict 
ourselves to the lowest order processes qq — > 7* — > e + e~, i.e., the electromagnetic 
annihilation of quarks and antiquarks. Contributions from QCD processes like qq — > 



ge + e~ and qg — * ge + e~ have been considered [24| and appear to be not very significant 
(typically they increase the yield by a factor less than 2.5). The rate of dilepton 
production per space-time volume can be obtained within a kinetic theory approach as 
dN e+e - f d 3 p 1 d 3 p 2 d 3 p + d 3 p_ 

h{X,Pl) J2{X,P2) 



d 4 x J 2E 1 (2tt) 3 2E 2 (2tt) 3 2£ + (2tt) 3 2E^(2tt) 

x\M(12 - e+e-)| 2 (2tt) 4 S^fa + p 2 - P + - P-), (3) 

where p a = (E a ,p a ) are the four momenta of incoming quarks (a = 1, 2) and outgo- 
ing leptons (a = e + ,e~). The quantity |.M(12 — > e + e~)| 2 denotes the square of the 
amplitude of the process, summed over initial and final state particles. The quantities 
fa{x,p a ) stand for the distribution functions of quarks and antiquarks. Integrating 



eq. (3) over the momenta of the leptons one gets a widely used expression [^5| for the 
dilepton rate per four- momentum of the pair, Q = p + + p_, 



dN e + e - r d 3 pi d 3 p 2 



h(x, Pl ) f 2 (x,p 2 )va(M 2 ) S (4 \ Pl +p 2 - Q). (4) 



d*xd*Q J (27r) 3 (27r) 3 
The pair's four momentum can be expressed as Q = (MichY, Q±, M±shY) with M± 



M 2 + Q\ as transverse mass, M as invariant mass, and Y as rapidity of the pair. 
The relative quark velocity reads v = M/(2EiE 2 ), while the total cross section is 
a(M 2 ) = |^12Fg with quark charge factor F q = J2 e l — § f° r u ;d quarks, and a is 
the fine structure constant. 

For the calculation of dilepton spectra in the region of large invariant mass one can 
utilize the Boltzmann approximation for the incoming partons 

f a (x,p) = X a (x) exp j P " T ^y I , ( 5 ) 

where u is the four- velocity of the medium. The above parton distributions are assumed 
to be in thermal equilibrium in momentum space, but not necessarily in chemical 
equilibrium. Only for chemical equilibrium the fugacities A a become unity. 



B. M± scaling 



The space-time evolution of the produced deconfmed matter is governed in our 
approach by the longitudinal scaling-invariant expansion accompanied by quark and 
gluon chemical equilibration processes |]2^j26| . Within such a model the space-time 
volume reads d 4 x = TrR^Tdrdfj with r as proper time and flow rapidity fj. Integrating 
the rate eq. (4) over the parton momenta and space-time evolution results in the 
dilepton spectrum 

«UW _«>Rl F J dTTK (Mi_\ 



dM\dQ\dY 4vr 2 q J \T{t 

The physical information encoded in the spectrum eq. (|6|) can be inferred from the 
following approximation. For large values of M± the main contribution to the time 
integral stems from early, hot stages at r ~ 7$. Under the assumption that the equation 
of state behaves like p oc e (here p and e stand for the pressure and energy density, 
respectively) one can derive from energy-momentum conservation the relation T = 
Tiiri/rY^Fir), F = [(d g Xf + d q Xf) / '(d g X 9 (r) + d 9 A 9 (r))] 1/4 with d q , g as degeneracy 
factors of quarks and gluons. Since near n one can approximate F ~ 1 and X q ~ Af, 
one arrives at 



dN e+e - 3a 2 R 2 A F q (nXfTfV u (M ± 



H hjr (7) 



dY dM\ dQ 2 ± 4tt 2 V M l / ^ T * 

with H(x) = x 3 (8 + x 2 ) K 3 (x) and K n as modified Bessel function of nth order. For 
the possible initial conditions at LHC (see below) eq. (^) represents a fairly good ap- 
proximation, also for small values of Xf 9 . Eq. ([?[) has the structure /i(T,Af) f 2 (M±/Ti), 
therefore one can infer from it the value of Tj by measuring the transverse rate at 
two distinct values of M±. Afterwards, the combination r^A^ can be extracted. If one 
could constrain by other means the initial time of the thermalized era, Tj, then even 
eq. (0) allows to estimate the initial fugacity Xf. Due to particle production processes 
the comoving entropy density does not longer serve as a link between initial and final 
states in the evolution dynamics. 

Our choice of initial conditions for produced deconfmed matter is based on the 



estimates of refs. [|T|,^7]] f° r ^ ne mm ij e t plasma with a suitable parametrization of the 



soft component which are similar to the self-screened parton cascade model [gc 
We take as main set of parameters for the initial temperature Tj = 550 (1000) MeV, 
for gluon fugacity Af = 0.5, and for light quark fugacity Xf = |Af of the parton plasma 
formed at RHIC (LHC) at initial time Tj = 0.2 fm/c. To achieve an upper "optimistic" 



limit of the thermal dilepton yield we consider below also higher values of the initial 
quark and gluon fugacities (Af = 0.5, 1, and Xf = 1) which lead to different final states. 
For the sake of definiteness we assume full saturation at confinement temperature T c = 



170 MeV and a quadratic time dependence of X q,9 (r) Hl5"l , [24 



To get realistic dilepton spectra which can be measured in relativistic heavy-ion 
collisions one should take into account the experimental conditions related to the ac- 
ceptance of detectors. To do so we employ Monte Carlo simulations with proper dis- 
tributions of created leptons with respect to individual rapidities y± and transverse 
momenta p±±. In case of the thermal emission of leptons from deconfined matter 
the needed distribution can be obtained directly from eq. (3). Due to the energy- 
momentum conservation the exponential function, originating from the Boltzmann 
distribution, can be rewritten as exp{(pi + p 2 ) ■ u/T} = exp{(j» + +p_) • u/T}. Since 
also |.M(12 — > e + e _ )| 2 = \A4(e + e~ — > 12) | 2 the integration in eq. (3) over the parton 
momenta yields 

dN e+e - f d?p + d 3 p_ j (p + +p_) -u\ + _ 

exp<^ }Fa{e + e 12), (8) 



d 4 x J 2£+(27r) 3 2£_(27r) 3 * \ T 

where cr(e + e~ — > 12) = <r(12 — > e + e~) = a(M 2 ) and the flux factor becomes F = 2M 2 . 
Then the space-time integrated distribution of thermal leptons reads 

dNe+e- o?R\ r ( M ± \ 2 

with dilepton transverse mass M]_ = p\ + + p\_ + 2pj_ + pj___ch(?/ + — y_). To get the 
invariant mass spectra we also use the kinematical relationship M 2 = 2pj_ + p_ s __[ch(y + — 
yJ) — cos(0 + — 0-)], where 0± denote the azimuthal angles of the leptons in the 
transverse plane. 

As pointed out above, due to the dominant contribution of correlated open charm 
and bottom decays via D(B)D(B) — > e + Xe~X into the high invariant mass dilep- 
ton spectra it is very unlikely that the thermal dilepton production from quark fu- 
sion can be observed by using the distribution dN e + e -/dM 2 alone even when impos- 
ing tricky cuts. In this situation the more detailed information of the full spectrum 
dN e + e - j dM\ dQ\ dY is expected to allow for better chances to extract the wanted 
thermal signal. 

As known since some time [|4],[^,[3(J and as directly seen in eq. (6) the dilepton 
yield from an equilibrium quark-gluon plasma depends only on the transverse mass, i.e., 
dN e+e -/dM 2 ± dQ\dY oc f(M ± ); that is, it scales with M ± . The main conditions for 
this famous M± scaling are p9f (i) local thermalization of the source of the dileptons, 



(ii) no scale other then the temperature, and (iii) predominant boost-invariant flow of 
matter. The most convincing way to see the validity of the Mj_ scaling is to consider the 
Q± dependence of the dilepton spectrum at fixed value of M±. In case of the boost- 
invariantly expanding quark-gluon plasma one gets for dN e + e - /dM]_\M ± =& x dQ\ dY 
indeed a straight line between Q± = and M±. A detailed analysis of the transverse 
momentum dependence of the dilepton spectra is performed in ref. [23], where the 
possible origins of a M± scaling violation related to transverse flow and contributions 
from the hadron gas are studied. The transverse expansion of deconfined matter affects 
hardly the M± scaling, while for hadron matter the life time is strongly reduced, so that 



its contribution to the total spectrum is moderate. In fact, our previous studies p4 



show that the conditions for the Mj_ scaling are satisfied for radiation from thermalized 
matter expected to be created in central collisions at RHIC and LHC. At the same time 
the dileptons from correlated open charm and bottom decay should strongly violate 
the Mi scaling if the cc and bb pairs are dominantly produced by hard initial parton 
collisions. The qualitative argument is as follows: Similar to the Drell-Yan process 
the hard charm and bottom pairs are created mainly back-to-back (unless hard gluon 
radiation spoils this correlation), i.e., the pair momentum is preferential ~ 0, 

and consequently also the Q± spectrum at fixed is strongly peaked at 0. The 

subsequent semileptonic decays cause a partial randomization of the lepton directions 
but, due to the Lorentz boost effect, both the leptons still are concentrated at small 
values of Q±. As a result one can expect that in the large Q± region the thermal signal 
sticks out, while the semileptonic decay products are sharply dropping. 



C. Monte Carlo simulations 



To quantify this effect we perform Monte Carlo calculations generating thermal di- 
electrons in accordance with the distribution eq. (9). As above we take into account the 
acceptances of PHENIX at RHIC and ALICE at LHC. First we consider the simplified 
case without any momentum degradation of heavy quarks, i.e., energy loss effects are 
not included and the 5 function fragmentation scheme is employed. This is obviously 
the most favorable case for the relative contribution from decay dileptons to the full 
dilepton yield. 

In order to get some threshold effect for the decay dileptons in dependence of 
at fixed Mi one needs to select leptons with individual transverse momenta greater 
than the maximum electron energy in the rest system of the decaying open charm 



or bottom meson. Typically such a threshold value is given by 1 (2) GeV for charm 
(bottom). The results of our calculations of the Q± spectra for RHIC conditions are 
displayed in fig. 3a. We select here dileptons with transverse mass within the interval 
2 GeV < Mi < 3 GeV. The value of the minimum transverse momentum of leptons is 
chosen as p™ m = 1 GeV. The detector acceptance destroys the ideal M± scaling, i.e., 
the thermal yield gets some structure reflecting the p± cut and the geometry. One can 
see the threshold like behavior of the dilepton spectrum from charm decays: there is 
a rapid dropping of the yield at Q± > 1.3 GeV. On the contrary the thermal dilepton 
spectrum sticks out up to the kinematical boundary Q± = M±. Unfortunately such 
a value of p™ in = 1 GeV is not large enough to cause a similar suppression effect for 
leptons from bottom decays. Nevertheless as pointed out in ref. J7[ and as seen in fig. 3a 
the contribution from bottom decays is suppressed (this is due to the narrow rapidity 
gap covered by PHENIX) but it is clearly above the thermal yield for the initial quark 
fugacity X q = 0.1. To estimate the chances to see a thermal dilepton signal over the 
background from bottom decays we perform variations of the initial quark fugacity. 
Due to the factor X 2 in the thermal rate eq. (9), even a comparatively small increase 
of the initial quark phase space saturation can result in a noticeable enhancement of 
thermal dileptons compared to the decay background. This can be seen in fig. 3a for 
the initial quark fugacities X q = 0.5 and 1. For the latter value the thermal signal is 
above the bottom decay contribution. 

For LHC conditions one can obtain the threshold like effect in the Q± spectrum 
even for the dielectrons from the massive bottom decays. In fig. 3b we show our results 
for dileptons within the transverse mass interval 5.25 < M± < 5.75 GeV. The electron 
transverse momentum cut p™ m = 2 GeV causes a sufficiently strong threshold effect 
both for charm and bottom, while the dielectrons from thermalized matter exhibit 
the approximate M± scaling even for such a value of p™ in . It gives therefore the 
opportunity to measure the thermal dilepton signal on the basis of the M± scaling 
at LHC conditions. Actually such measurements could be performed for lepton pairs 
with high enough transverse momentum but smaller invariant mass. In particular, as 
follows from fig. 3b, for a fixed value of M± ~ 5.5 GeV the proper region for transverse 
pair momenta is about 4 GeV < Q± < 5.5 GeV. This window increases considerably 
with increasing initial quark fugacity A* of thermalized parton matter. 

We perform also calculations of the Q± spectra at fixed M± for decay dileptons 



including the Peterson fragmentation function |[31|| , where the function 



/(*) 



z{l-- 

z 



(10) 



(e = 0.06 (0.02) for charm (bottom) [p2jl ) prescribes the probability to change the 
spatial momentum of a heavy quark from p to zp during hadronization. Independently 
we also study the influence of energy loss of heavy quarks in the deconfmed medium. 
Both these effects cause a suppression of the dilepton yield from open charm and 
bottom, in particular in the small-Qj_ (large-M) region. Since the suppression effect 
due to Peterson's fragmentation is quite small (typically within factor 2) we do not 
display it separately. On the contrary, the degradation of the heavy quark momenta 
due to energy loss in the deconfmed medium causes a considerably suppression of the 
contribution of decay leptons in the Q± spectrum at fixed value of M±. In fig. 4 we 
show the results of our calculations for LHC conditions with the same kinematical cuts 
as in fig. 3b. As seen in this plot, both charm and bottom are considerably suppressed 
due to energy loss effects, in particular in the region of smaller values of Q±. These 
calculations exploit the model I. In case of model II the suppression effect is much 
stronger. With respect to the wanted thermal dilepton signal the suppression of the 
Q± spectra from bottom decays appears most important. The proper window to 
extract thermal dileptons becomes greater, e.g., 3 GeV < Q± < 5.5 GeV for initial 
quark fugacity X q = 0.5. Therefore, the energy loss effect obviously helps to measure 
thermal dielectrons with ALICE at LHC on the basis of the M± scaling behavior. 

To clarify whether the proposed signal is experimentally feasible we estimate the 
background from 7r° decays within the same kinematical gates. For this background 
study we also perform Monte Carlo simulations of tt° — > e + e~7 decays utilizing the 
distribution dN n o/dydp± = (dN^o/dy) exp{px/T n } with T n = 200 MeV in accordance 



with present SPS results for central Pb + Pb collisions at 158 AGeV beam energies [33 



The slope parameter T n is a convolution of local freeze-out temperature and transverse 
flow. The latter one is expected to increase from SPS to RHIC to LHC due to higher 
rapidity densities, and it might cause a slightly higher value of T n . We choose here 
the fictitious values dN^o/dy = 2500 (500) for LHC (RHIC). Each vr° is three-body 
decayed in its rest frame with the branching ratio of 1.2% |Z| and then the resulting 
e + e~ pair is boosted and filtered. We find the individual 7r° decay dileptons with 5.25 
GeV < Mj_ < 5.75 GeV three orders of magnitude below the thermal dilepton signal 
and only at Q± > 5.2 GeV. The flattening of the pion spectrum at large values of 
p_L, due the Cronin effect and hard scatterings, is not covered by our parametrization 
of the p± spectrum; it will enhance the pion decay contribution. Otherwise some 



quenching |T2j at large values of p± can counteract and balance such effects. The finite 
experimental momentum resolution will cause a decay contribution also at somewhat 
smaller values of Q±. However, these effects will hardly close the window for the 
thermal signal. The next-to-leading order Drell-Yan contribution is also estimated to be 
strongly suppressed at larger values of Q± (say, for Q± > 1 GeV) because of the chosen 
large p± cut. The most prominent danger however comes from the combinatorical 
background from uncorrelated hadron decays: Only after experimentally removing 
this background together with the background from the conversion dileptons, emerging 
from energetic photons in the detector material, one has a chance to see the interesting 
signal. 

Fully analog considerations also apply at RHIC. However, selecting electron pairs 
with such large values of M± pa 5.5 GeV and choosing p™ in = 2 GeV results in pro- 
hibitory small counting rates which seem to make such a strategy not feasible. For 
the above choice of M± pa 2.5 GeV and choosing p™ m = 1 GeV the dileptons from 
7T° decays produce here a contribution which is three orders of magnitude above the 
thermal signal, however only in the interval Q± > 1.95 GeV. Therefore, one remains 
in a situation as depicted in fig. 3a and discussed above. 

IV. CONCLUSIONS 

In summary we analyze the chances of observing a thermal dilepton signal in high- 
energy heavy-ion collisions. Exploiting currently estimated energy loss effects of heavy 
quarks passing through deconfined matter we find that the invariant mass spectra of 
dileptons, as measurable in the detector facilities PHENIX and ALICE, offer hardly a 
chance to find thermal dileptons with M > 2 GeV. Even if in optimistic energy loss 
scenarios the dileptons from correlated charm and bottom decays compete with the 
Drell-Yan background, the thermal signal is significantly smaller. Since the energy loss 
is a poorly quantified effect we stress the need to make attempts to verify it, for in- 
stance via the single electron p± spectra which should reflect the transverse momentum 
degradation of parent charm and bottom mesons. 

The double differential rate dN e + e - / dM\dQ\ with M± in a narrow interval and 
with a suitable p™ in cut on the individual leptons, however, seems to allow for a window 
at large values of the pair transverse momentum where the thermal yield shines out. At 
RHIC this is probably not feasible, while at LHC the estimated high initial temperature 
of the parton matter causes a stronger thermal signal. The basic issue is here whether 



the combinatorical background of many other strong sources can be subtracted with 
sufficient accuracy. 
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FIG. 1. The transverse momentum spectra of single electrons from D and B meson decays 
(a) and their sum (b) at RHIC energies within the PHENIX acceptance. Displayed are the 
spectra before ("initial") and after energy loss within model II. 
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FIG. 2. Dielectron spectra from correlated charm and bottom decays as a function of the 
invariant mass; long-dashed (solid) lines are for results without (with) energy loss according 
to model II for RHIC with PHENIX acceptance (a) and model I for LHC with ALICE 
acceptance but without high-p± cut (b). The dashed lines depict the thermal dielectron 
spectrum. For the strength of the Drell-Yan contribution cf. [9,10] 
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FIG. 3. The transverse momentum spectra of dielectrons at constrained transverse masses: 
(a) RHIC with PHENIX acceptance, p± > 1 GeV and 2 GeV < M± < 3 GeV, (b) LHC 
with ALICE acceptance, p±_ > 2 GeV and 5.25 GeV < M± < 5.75 GeV. The dash-dotted 
(short-dashed) lines are for correlated charm (bottom) decays. The full, long-dashed and 
gray-dotted lines depict the thermal yield for various initial quark and gluon fugacities. 




FIG. 4. The same as in fig. 3b but with energy loss effects according to model I. 



